Search for extraterrestrial point sources of high energy neutrinos with AMANDA-II 

using data collected in 2000-2002 
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^ . The results of a search for point sources of high energy neutrinos in the northern hemisphere 

using data coUected by AMANDA-II in the years 2000, 2001 and 2002 are presented. In particular, 
a comparison with the single-year result previously published shows that the sensitivity was improved 
by a factor of 2.2. The muon neutrino flux upper limits on selected candidate sources, corresponding 
to an neutrino energy spectrum, are included. Sky grids were used to search for possible 

excesses above the background of cosmic ray induced atmospheric neutrinos. This search reveals no 
statistically significant excess for the three years considered. 

PACS numbers: 95.85.Ry,95.55.Vj,96.40.Tv,98.54.-h 



The detection of high energy cosmic rays raises funda- 
mental questions about their generation and the mecha- 
nisms responsible for such energies. The origin of cosmic 
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rays above the "knee" (10^^ eV) still remains uncertain. 
Nevertheless there is evidence that below such energies 
they are generated through acceleration mechanisms in 
expanding supernova remnant shocks H and in micro- 
quasars 0i 0; Bjj although wc cannot exclude the pos- 
sibility of extragalactic sources at these energies. The 
interaction of accelerated protons with ambient matter 
or radiation leads to pion production and, consequently, 
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to neutrinos and gamma rays following a power-law en- 
ergy spectrum. High energy gamma rays are affected by 
absorption during propagation and may also be produced 
by inverse Compton scattering of shock-accelerated elec- 
trons. Therefore detection of gamma rays alone is in- 
sufficient evidence for hadronic acceleration. Neutrinos 
can provide a link to increased understanding of high 
energy cosmic rays, although most predictions of high 
energy extraterrestrial neutrino fluxes conservatively re- 
quire kilometer-scale detectors 0, S • 

AMANDA-11 (Si operates at the geographic south pole. 
It is composed of 677 optical modules (OMs) - photomul- 
tiplier tubes encased in glass pressure vessels - spaced 
along 19 vertical cables (strings) arranged in concentric 
circles. The instrument spans a geometrical volume of 
clear glacial ice between depths of 1500 and 2000 m, 
with a diameter of 200 m. The AMANDA-II neutrino 
telescope has been in operation since 2000. In this let- 
ter, we follow up on a previously published search for 
high-energy neutrino point sources from the data sam- 
ple collected in 2000 3, using data from the three years, 
2000 to 2002. The sensitivity for the detection of point 
sources has significantly improved in AMANDA, com- 
pared to 1997 [13 and 1999 results 0|, due to both 
detector performance and analyses technique improve- 
ments. 

A high-energy muon neutrino interacting with the ice 
or bedrock in the vicinity of the detector produces a high- 
energy muon propagating a few kilometers when above 
1 TeV. At these energies the mean angular offset be- 
tween the muon track and incident neutrino is less than 
0.8°. The muon track is reconstructed using the detec- 
tion of Cherenkov light emitted as it propagates through 
the array of OMs, and the likelihood of arrival time of 
detected photons at each OM location. The resulting 
zenith-dependent median pointing resolution varies be- 
tween 1.5° and 2.5° [12. 

Muons are also produced by the interactions of cosmic 
rays in Earth's atmosphere. These atmospheric muons 
dominate the AMANDA-II trigger rate so the search 
for neutrino-induced muon tracks is only conducted in 
the northern hemisphere, using Earth as an atmospheric 
muon filter. A second source of background is repre- 
sented by atmospheric muon tracks reconstructed as up- 
going. These events can be rejected using track quality 
criteria. The most important source of background is the 
residual up-going flux of neutrinos produced in the atmo- 
sphere by the impact of cosmic rays. These atmospheric 
neutrinos also serve as a verification of the detection prin- 
ciple and demonstrates AMANDA's capability as a neu- 
trino detector 0, 0| . The search for possible extrater- 
restrial neutrinos begins with a dataset dominated by the 
well- understood atmospheric neutrinos. This analysis se- 
lects a three- year sample of events with median energy 
of ~ 1.3 TeV and extending up to ~ 100 TeV. Extrater- 
restrial neutrinos are beleived to be distinguished by a 
harder energy spectrum, taken as proportional to E^'^ in 
this analysis. 



The exposure of the present analysis is three times 
higher than that of the previous analysis [3- A differ- 
ent search strategy is used, which includes an explicit 
high energy event selection to reduce the expected lower 
energy atmospheric neutrino background. 

I. DATA ANALYSIS 

The data used for this analysis were collected between 
the months of February and November in the years 2000, 
2001 and 2002 (see Table I). 



year 


livetime (days) 


triggers 


2000 


197 


1.34 X 10'-' 


2001 


194 


2.04 X lO'* 


2002 


216 


2.17 X lO'' 



TABLE I: The experimental livetime and number of triggered 
events for each year used in this analysis. The triggered events 
may vary in different years mostly due to different cleaning 
procedures, which are mainly affected by the number of stable 
OMs during the specific year. 

The experimental sample used in this analysis corre- 
sponds to a total of 607 days of livetime and contains al- 
most 5.6 billion triggers. Starting from 2002, a first level 
filter is performed at the South Pole during data taking. 
The reduced amount of data is transferred via satellite to 
the northern hemisphere for analysis. After the applica- 
tion of an iterative maximum-likelihood reconstruction 
algorithm and the selection of tracks that are likely to 
be upgoing (l^, about 0.45 million events with recon- 
structed declination S > —10° remain. Since AMANDA- 
II is located at the south pole, S = 0° corresponds to hori- 
zontal and 6 = 90° to vertical up-going directions. These 
events, containing mostly mis-reconstructed atmospheric 
muons and a contribution of atmospheric neutrinos, were 
used as an experimental background for selection opti- 
mization. 

To avoid biasing the event selection the data were 
scrambled by randomizing the reconstructed right as- 
cension (a) of each event. The optimization procedure 
makes use of three observables: the number of hit OMs 
for each event (nch), the reconstructed track length in the 
array and the likelihood ratio between the muon track re- 
construction and a muon reconstruction constrained by 
using an atmospheric muon prior ^ full simula- 

tion chain, including neutrino absorption in the Earth, 
neutral current regeneration, muon propagation and de- 
tector response for the given data taking periods, is used 
to simulate point sources of muon neutrinos and anti- 
neutrinos 9] . Events are simulated at the center of each 
5° band of declination (S), according to an energy 
spectrum. The final cuts on these observables and the 
optimum size of each circular search bin were indepen- 
dently determined for each declination band in order to 
have the strongest constraint on the signal hypothesis. 
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This corresponds to the best sensitivity, i.e. the average 
flux upper hmit obtained in an ensemble of identical ex- 
periments assuming no signal The true directional 
information was then restored for the calculation of the 
limits. 

The upper limits of this analysis were calculated using 
the background nt measured using the events off-source 
in the corresponding declination band, and the expected 
number of events, Ug, from a simulated point source of 
known flux ^{E): ^umit{E) = ^{E) x ^90 (riob^, nb)/^^. 
Here Uots is the number of observed events in the given 
source bin, and /xgo is the upper limit on the number 
of events following the unified ordering prescription of 
Feldman and Cousins 0| . The three years were analyzed 
both separately and as combined data samples. 
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III. RESULTS 

Figure ^ shows the calculated sensitivity versus dec- 
lination for energies above 10 GeV. The event selection 
used produces a sensitivity which is fairly constant over 
all dechnations. For 0° < 5 < 5° {0 < sm{5) < 0.09) 
the background contamination is 4 times higher than for 
(5 > 5°, and the sensitivity is poorer. For 6 > 80° (sin((5) 
> 0.98), on the other hand, the solid angle aperture is 
small and the background evaluation is affected by higher 
relative statistical fluctuations. The sensitivity improves 
with the detector exposure and for the three-year sample 
it is ~ 2.2 times better than for a single year. This im- 
provement is better than would be expected from longer 
exposure alone in the presence of background, due to im- 
proved analysis techniques. 



year 








2000f9] 


601 


676 


133 


2000 


306 


296 


111 


2001 


347 


364 


115 


2002 


429 


429 


131 


00-02 


646 


635 


297 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin(5) 



TABLE IL The number of observed events with 5 > 5° after 
cut optimization, for each year and the combined three-year 
sample. The numbers relative to reference are compat- 
ible with a normalization factor of ~ 0.86, for the atmo- 
spheric neutrino simulation, as quoted in the above refer- 
ence. The numbers Up of the predicted atmospheric and sig- 

nal neutrino events (with signal energy spectrum of ~ 



FIG. 1: Sensitivities on the integrated flux above i?^ = 
10 GeV as a function of declination and for an energy 
spectrum. The sensitivities for the year 2000, 2001 and 2002 
are compatible with each other, and shown along with the one 
for 2000-01 and for the 2000-02 three-year sample. 



II. CALIBRATION AND SYSTEMATIC 
UNCERTAINTY 

Using the three-year experimental sample the abso- 
lute normalization of the detector simulation, with re- 
spect to atmospheric neutrino flux, was determined to be 
1.03 ± 0.02. This normalization factor is different from 
the value in reference but consistent with the theo- 
retical uncertainty of 25% The optical properties 
of the fiducial ice, needed for the detector simulation, 
are determined using down-going muon data and in-situ 
calibration lasers [13 • The overall experimental system- 
atic uncertainty in the acceptance was evaluated to be 
~ 30% [23 ■ The absolute pointing accuracy, determined 
with coincident events between the SPASE air shower ar- 
ray lU and AMANDA-II, is better than 1°, i.e. smaller 
than the angular resolution. These systematic uncertain- 
ties do not affect significantly the limit calculations (^] . 



The final three-year sample consists of 646 upward 
{6 > 5°) reconstructed muons (see Table II). The pre- 
dicted number of atmospheric neutrinos is 635. In the 
year 2000 alone, the number of selected events is 306, 
compared with the 601 (699 for (5 > 0°) in reference j^. 
The difference between the two samples is due to the 
different choice of observables used for the selection op- 
timization. In particular the use of the nch observable, 
which is correlated to the energy released by the muon in 
the array and, ultimately, to the neutrino energy, selects 
~ 26% higher median energies than those in [3 (from 
~ 700 GeV to ~ 1 TeV for a single year). This selection 
is obtained at the price of removing a significant frac- 
tion of atmospheric neutrino events: for instance only 
221 events in [3 would survive the new selection, 94% of 
which (i.e. 207) are also found in the new sample from 
the year 2000. 

As shown in Table II the number of events in the fi- 
nal sample (nofcs) docs not sum up with experimental 
exposure. Since the signal hypothesis predicts a higher 
event intensity at high energy than the atmospheric neu- 
trino background, a longer exposure allows a stronger 
constraint on a given model by requiring a stronger en- 
ergy cut (the median energy increases from ~ 1 TeV for a 
single year to '--^ 1.3 TeV for the three-year sample), which 



rejects more background events and results in stricter 
limits. Consequently the three-year sample contains ~ 
40% fewer observed events than the sum of single years, 
but only ~ 17% of the high energy neutrino signal events 
are lost. 
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FIG. 2: Muon neutrino effective area as a function of the neu- 
trino energy at different declinations. The effect of neutrino 
absorption in the Earth is responsible for the effective area 
decrease at high energies and declinations. 

The detector performance is assessed by the neutrino 
effective area A'^jj:{E,y,S), which contains the neutrino 
interaction probability, muon propagation, detector re- 
sponse and the analysis selection. It is defined by the 
relation between the differential neutrino flux .i'^'i, and 
the predicted number of neutrino events np{iy) 
the equation 
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Figure 13 shows the muon neutrino effective area as a 
function of neutrino energy for the three-year optimized 
selection. The curves are shown for different declina- 
tions. Above 1 PeV neutrinos begin to be absorbed by 
the Earth, except for the events that enter AMANDA-II 
horizontally. 

A binned search for excesses in the b° < 5 < 85° 
region was performed on the three-year event sample. 
The search grid contains 290 rectangular bins with 
declination-dependent width ranging from 5.6° to 8.8°, 
based on the optimized search bin diameter. The grid is 
shifted 4 times in 5 and a to fully cover boundaries be- 
tween the bins of the original configuration. A higher 
number of grid shifts showed no improvement in the 
average maximum statistical significances on simulated 
Poisson-fluctuated signal with intensities comparable to 
the background. The probability distribution for back- 
ground fluctuations in the ensemble of bins was evaluated 
by using 20,000 experimental samples with scrambled a 
and calculating the highest value of the maximum statis- 
tical fluctuation significance over the entire sky. 



TABLE III; 90% CL upper limits on candidate sources. Re- 
sults from the present analysis are reported for a comparison 
with the limits from 1^. Limits are for the assumed E~'^ 
spectral shape, integrated above Ev ~ 10 GeV, and in units 
of lO-^cm-^s-i 



The bin with the most statistically significant excess 
from the three-year experimental sample is at about 
a = 22h and S = 21°, with 10 observed events in the 
search bin on a background of 2.38 events, estimated from 
the corresponding declination band. The observed excess 
has a statistical significance of 1.9x10"^ (3.73 a). The 
chance probability of such an excess, in the ensamble of 
bins, is 28%. 

Table III shows the 90% CL neutrino flux limits for 
northern hemisphere TeV blazars, selected GeV blazars, 
microquasars, magnetars and selected miscellaneous can- 
didates. The limits are compared with the values from 
0, they are compatible with the average flux upper limit, 
or sensitivity, of Figure ^ and the deviations from it are 
due to statistical fluctuations in the observed sample. 

Figure 13 shows the 90% CL neutrino flux upper limits 
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FIG. 3: 2000-02 upper limits (90% CL) on the neutrino 
flux integrated above 10 GeV in equatorial coordinates for 
5 > 5°. Limits (scale on right axis) are given in units of 
xlO^* cm~^s~^ for the assumed spectrum. Systematic 
uncertainties are not included. The cross symbols represent 
the observed events. 

in equatorial coordinates. The limits are calculated by 
scanning the sky and counting the events within the op- 
timized search bins at the given declination. The highest 
upper limit in the Figure corresponds to the previously 
discussed statistically significant bin. Other high limit 
spots visible in the figure have statistical significances 
smaller than 3.4 cr. 

We analyzed the 2000-02 data sample collected by the 
AMANDA-II detector to search for point sources of high 



energy neutrinos. We performed both a non-targeted 
binned search and a targeted search focussing on known 
objects that are potential high energy neutrino emitters 
(as in reference Q). We found no evidence of a signifi- 
cant flux excess above the background. A km-scale ex- 
periment, such as IceCube |22|, will be able to increase 
the detection sensitivity by at least a factor of 30 in the 
same time scale. 
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